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ABSTRACT 



A computer program to design geosynchronous spacecraft has been developed. The program 
consists of four separate but interrelated executable computer programs. The programs are compiled 
to run on an dos based personnel computer. The source computer code is written in DoD mandated 
Ada programming language. 

The thesis presents the design technique and design equations used in the program. Detailed 
analysis is performed in the following areas for both dual-spin and three axis stabilized spacecraft 
configurations: 

• Mass Propellent Budget and Mass Summary 

• Battery Cell and Solar Cell Requirements for a Payload Power Requirement 

• Passive Thermal Control Requirements 

Thesis includes a users manual Appendix A, and the source code for the computer programs 
as Appendix B. 
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I. INTRODUCTION 



A. PURPOSE 

The purpose of this thesis was to design an executable computer program capable of 
determining the numbers required to design specified aspects of a geostationary 
communications satellite, that could be run on a standard home personal computer 
(PC) The program is based on Professor Brij N. Agrawal’s book "Design of 
Geosynchronous Spacecraft". The book describes all the steps necessary to design a 
geosynchronous communications satellite. Because of time limitations, it was decided 
that the following areas of a Geostationary Communications Satellite would be 
completed on this portion of the project: 

• Mass and Propellent Budget 

• Electric Power 

• Thermal Control 

Both Three-Axis and Dual Spin Stabilized satellite configurations are included in the 
program and provide design reports that can be printed at the user’s discretion. 
Appendix A is a short user’s manual of the computer program and Appendix B 
contains the different executable computer programs’ source code. 
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B. BACKGROUND / CONCEPT 



The program came about through discussions with Professor Brij N. Agrawal 
and space system students enrolled in his satellite design classes. The basis for the 
computer program is to aid the user in quickly determining useful and accurate 
numbers for a geostationary spacecraft design. Although not a tutorial, the program 
sequentially walks the user through the necessary steps to develop a Mass Propellent 
Budget, Photovoltaic Electric Power System, and Thermal Control System for either a 
three-axis or dual-spin stabilized geostationary communications satellite. The four 
executable programs will, with a small amount of preparation, allow the average user 
to develop the described areas of a geosynchronous satellite design in minutes. 

Where as by using a hand calculator, the same variable numbers and final results 
would take days or weeks to arrive at the same results. 

C. OPERATIONAL DESIGN PARAMETERS 

All parameters for the executable program used to design a spacecraft in this 
thesis are based on an orbiting geostationary communications satellite. Some velocity 
components for propellent usage and therefore mass budget are related to the transfer 
and parking orbits as well as orbit injection angles and delta velocities necessary to 
achieve and maintain geostationary orbits. The equations or components used are 
specifically applicable to a geostationary satellite and might be useful for others as 
well. The two most common satellite design types, three axis stabilized and dual-spin 
stabilized, are used as a basis for program engineering and concept. 
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1. Dual Spin Stabilization 

The dual spin configuration used for executable program design is a 
cylindrical shell with spacecraft power provided by solar cells surface mounted around 
the outer cylindrical area. The interior of the shell of the cylinder houses spacecraft 
electronics and propulsion devices. 

2. Three Axis Stabilization 

Three axis stabilized spacecraft will be based on a rectangular box shape. 
Power will be provided by sun tracking flat panel array(s). During parking orbit and 
operational orbit insertion, spin stabilization will be used. 
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II. MASS SUMMARY AND PROPELLENT BUDGET 



This chapter will step through the process of determining a equipment mass 
summary and a propellent weight budget. 

A. TRANSFER ORBIT 

For a satellite to be useful it first must achieve its operational orbit. The 
spacecraft will progress through a series of different orbits before achieving its "final" 
working orbit. The sequence begins with the launch of the spacecraft from any of a 
variety of locations depending on payload type, mass, desired orbit inclination, including 
mission or launch window constraints. The satellite is released from the launch vehicle 
into a Low Earth Orbit (LEO) parking orbit where system checks are performed. The 
next phase is to put the satellite into an elliptical transfer orbit in preparation for final 
insertion into operational orbit. 

There are a variety of "transfer" orbits that can be used for final insertion into 
operational orbit. The most common in-plane and fuel efficient transfer orbit for a 
geostationary spacecraft is a Hohman Transfer Orbit (HTO). The programs in Appendix 
B use HTO for velocity change requirements for the transfer orbit. After release from 
the launch vehicle the spacecraft enters a parking orbit where satellite systems are 
checked out in preparation for insertion into the final operational orbit. 

The first phase of a Hohman transfer is to fire the perigee kick motor (PKM). 
This places the spacecraft in an elliptical transfer orbit. After four or five transfer orbits 
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the satellite completes system checkout, attitude determinations are finalized, and the 
apogee kick motor (AKM) fires. Geostationary orbit (GSO) is achieved by inserting 
the satellite into a circular orbit at a radius of approximately 42,160 kilometers and at 
a desired inclination (i) of zero degrees (equatorial orbit) for our design purposes. The 
spacecraft will then go through a series of re-orientations to finalize its operational orbit 
and mission altitude. (Wertz, 1991, p. 130 ) 

B. VELOCITY 

Velocity determination and certain delta velocities required to achieve the 
operational orbit are vital to building propellent summaries. Other factors taken into 
consideration are the type of fuel used (relating to I sp ) and the efficiencies of the separate 
orbit insertion and orbit maintenance activities. The efficiencies of the apogee kick 
motor AKM, PKM, station keeping, and de-orbit functions of the spacecraft must be 
known so that sufficient propellent, and propellent weight, is incorporated into spacecraft 
design. Margins for the propellent budget are also assimilated into the propellent budget 
as a safety factor even though the delta velocity equations provide very accurate figures 
with little error. All values determined are for a geostationary communication satellite. 

1. Geostationary Orbit 

In a perfect geostationary orbit the satellite will move around the earth’s 
equator synchronized with the rotation of the earth about the earth’s axis. The period 
(P) is equal to one sidereal day (23 hours, 56 minutes, 4.09 seconds), or from the rising 
of a star to the rising of the same star vice 24 hours or sunrise to sunrise. 
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The radius for a geostationary orbit is determined via Keplar’s Law using 
a P of 86, 164.09 seconds corresponding to one sidereal day and is given by 

a = ( ^ e * p2 ) 1/3 EQ 2.1 

4 *7I 2 

where n e = gravitational constant 398,601.2 km 3 /(kg/s 2 ) 
a = semi-major axis 
P 2 = period of the orbit in seconds 

The resulting orbital radius is 42,164.2 kilometers and is a pure circular orbit with zero 
eccentricity. 

2. Satellite Drift 

Any deviation from the above ideal parameters will cause perturbations 
in the orbit. The gravitational forces of the sun and moon also act on the satellite to 
cause drift. Table 2. 1 lists the average drift rate imposed on the satellite by the sun and 
moon. 



Table 2.1 SECULAR RATES FROM SUN AND MOON 





Effect of Moon 


Effect of Sun 


n 


-0.00076 


-0.00034 


co ddot 


*0.0 


*0.0 
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a. East - West Drift 

For a radius greater than (>) 42,164.2 kilometers the spacecraft’s 
orbital period is greater than 24 hours giving it the appearance of drifting westward 
with respect to the surface of the earth. Conversely, for an orbital radius of less than 
42,164.2 kilometers, the spacecraft’s orbital period is less than 24 hours, thus giving 
it the appearance of drifting eastward with respect to the surface of the earth. The 
incremental drift rate of the longitudinal drift with respect to an ideal geostationary 
orbit is: (Agrawal, pg 68, 1986) 

ft n ,- 3 = - 3 . 360°/ day _ -0 . 013°/ day eq 2.2 

2 a B 2 42,164. 2km 42,164.2 km 

The east - west motion for an inclined circular orbit is usually very small. This 
longitudinal drift for a small i and X is: 

X = - — *sin (2nC) EQ 2.3 

4 

So for small inclinations, latitude oscillation is predominate. This can be easily seen 
in Figure 2.1. The effects of the sun and moon can be seen in Table 2.1. For an 
elliptic equatorial orbit with inclination (i=0) the amplitude is 

&X=2*e EQ 2.4 

Therefore, east - west oscillation predominates. 
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Figure 2.1 Geosynchronous Orbit (i^O) 

b. North South Drift 

An inclined geostationary orbit with correct radius will tend to 
oscillate mainly in latitude with respect to the surface of the earth. The amount of 
oscillation depends on the inclination and the drift rate. The motion or path of the drift 
resembles a figure eight ( 8 ). See figure 2.1. (Agrawal, p. 88, 1986) 

c. Gravitational Effects of Sun and Moon 

The gravitational forces of the sun and moon cause secular 
variations in the orbital elements of a geostationary satellite. These variations have the 
largest effect on right ascension of the ascending node and the argument of perigee. For 
nearly circular orbits (one aspect of an ideal geostationary orbit) e = 0 and the 
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resulting error is on the order of e 2 . So the equations for the rates of change from the 
sun and moon are (Wertz, p. 125, 1991): 



• 00138* EQ 2.5 

fl sun = -0 • 00154* c °sU) EQ 2.6 



6 jnoon =-0 . 00169* (4 -5 *sin 2 (i) ) EQ 2.7 



(j £un =0. 00077* 



(4-5*sin 2 (i) ) 

n 



EQ 2.8 



where, i = inclination 

n = orbit revolutions per day 

For a geostationary orbit, the total longitudinal drift acceleration 
X is obtained by averaging the drift acceleration over the orbital period. Since the sun 
for our purposes is an inertially fixed reference point it contributes less than one degree 
of drift per day. So for a circular orbit the velocity remains constant. Therefore the 
main contributor to longitudinal perturbations is not the sun or moon but the earth’s 
slightly elliptical shape. This ellipticity caused by what is known as the equatorial 
bulge changes the earths gravitational force preventing it from acting purely radial. This 
phenomenon creates two stable longitudes called S 1 and S2 (75°E and 255°E) and two 
unstable longitudes US1 and US2 where the gravitational forces are radial and the 
longitudinal drift acceleration is zero. Drift will occur whenever the satellite is not at 
one of these points. For example, if a satellite longitude is between SI and US1, the 
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lateral force component is along the velocity and this results in a negative longitude 
acceleration as shown in EQ 2.9. So whenever is near US1 or US2 it will migrate 
towards the stable longitudes SI or S2 whichever is closest. (Agrawal, p. 83, 1986) 

i-0.0168.sin HX-X,' M ,)4SSI£22. EQ 2.9 

Taking into account only second order gravity effects the 
longitudinal drift acceleration is: 
where 

X = longitudinal drift acceleration in degrees \ day 2 
X = longitude of the satellite in degrees 
X stibl( .= stable longitude of 75°E or 255°E 

C. VELOCITY DETERMINATION 

Variables from the program strategy that will begin our design are based on a 
geostationary orbit. First we will assume that some launch vehicle has inserted our 
spacecraft into a desired parking orbit . For a final GSO, a good parking orbit might be 
around 6600 kilometers. A parking orbit is a LEO orbit where satellite system checks 
are performed prior to the other phases of insertion into the final operational orbit (GSO 
operational orbit — 42,160 kilometers). 

1. Orbit Insertion Velocities 

For the design calculations and for the executable program in appendix B 
the following equations are used for velocity determination. First, to find the orbit radius 
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the program requests the values for radius at apogee (r*) and then the value for radius 
at perigee (r P ) to find the semi-major axis (a). 



a = 



2 



EQ 2.10 



The equations for transfer orbit velocity at apogee (VJ and velocity 
transfer orbit velocity at perigee (V,,,) are: 



V, 



parking \ 



Pe 

a 



v = 

v t P ^ 



Z*Pe* r A 



V = y * ( f ) 

v ta v tp ' - ' 



EQ 2.11 



EQ 2.12 



EQ 2.13 



The in-plane velocity change required to transition from parking 
orbit to transfer orbit is: 



L.V r =V -V 

" transfer v tp Parking 



EQ 2.14 



and the mean velocity is given by: 




EQ 2.15 
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The transfer orbit period is: 



r t =2*7t*— 



i 



2 

2 




EQ 2.16 



And ultimately the desired synchronous orbit velocity is: 




EQ 2.17 



Delta Velocity for Apogee Motor Firing 



To find the velocity change and insertion angle required to 



transition from a parking orbit to the final geostationary equatorial orbit used in the 
program, the orientation of the plane of orbit or inclination must be changed. Figure 2.2 
depicts this change. The components of the velocity necessary to achieve equatorial 
geostationary orbit can be found using the law of cosines. The insertion angle using this 
law is : 



determine variable values. The HTO uses a two bum method for insertion into 
Geostationary orbit. This method is the most fuel efficient and is therefore highly 



a=arctan 



y ta *sin(i) 



EQ 2.18 



v e -^ ta *cos(i) 



A Hohman Transfer Orbit (HTO) is utilized by the program to 



12 



